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0^ \ We perform an O(a^) perturbative calculation of the equation of state of cold but 

dense QCD matter with two massless and one massive quark flavor, finding that 
perturbation theory converges reasonably well for quark chemical potentials above 
O ■ 1 GeV. Using a running coupling constant and strange quark mass, and allowing for 

^ . further non-perturbative effects, our results point to a narrow range where absolutely 

D I stable strange quark matter may exist. Absent stable strange quark matter, our 

findings suggest that quark matter in compact star cores becomes confined to hadrons 
only slightly above the density of atomic nuclei. Finally, we show that equations 

■ of state including quark matter lead to hybrid star masses up to M ~ 2M Q , in 
, agreement with current observations. For strange stars, we find maximal masses 

qq I of M ~ 2.75Mq and conclude that confirmed observations of compact stars with 

M > 2Mq would strongly favor the existence of stable strange quark matter. 
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I. INTRODUCTION 



The properties of cold nuclear matter at densities above that of atomic nuclei, in par- 
ticular its equation of state (EoS) and the location of the phase transition to deconfmed 
quark matter, remain poorly known to this day. The difficulty in performing first principles 
calculations in such systems can be traced back to the complicated non-linear and non- 
perturbative nature of Quantum Chromodynamics (QCD). These properties have precluded 
an analytic solution describing confinement, while non-perturbative numerical techniques, 
such as lattice QCD, are inapplicable at large baryon densities and small temperatures due 
to the so-called sign problem. This should be contrasted with the situation at small baryon 
density and large temperatures, where close to the deconfinement transition region lattice 
QCD has provided controlled results for the EoS as well as the nature of the transition [H, 0], 
while at temperatures much above the transition, the system is well described by analytic 
results from resummed perturbation theory 

Experimentally, the high temperature / low baryon density regime of QCD can be studied 
in relativistic heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) (7-l(| and 
in the future at the Large Hadron Collider (LHC) [ll|. Collisions at lower energy, e.g. at 
the Alternating Gradient Synchrotron (AGS) and Super Proton Synchrotron (SPS) [l2|, 



13] . as well as those planned at the Facility for Antiproton and Ion Research (FAIR) and 



RHIC [14|, [15|], study QCD matter at somewhat higher baryon density, and may give some 
insight into the EoS of cold nuclear matter. However, at truly low temperatures and supra- 
nuclear densities, QCD matter exists only in somewhat inconveniently located 'laboratories': 
Compact stars. 

In the cores of compact stars, nuclear matter is expected to reach densities several times 
that of atomic nuclei n sat ~ 0.16 fm~ 3 , so that astrophysical observations may be able to 
provide critical information about the EoS of strongly interacting matter in a regime inac- 
cessible to terrestrial experiments. Theoretically, the bulk properties of nuclear matter at 
or close to n sat have been studied using microscopic calculations [16] as well as phenomeno- 
logical mean- field theory [13]. While giving matching results for symmetric nuclear matter 
(equal number of protons and neutrons) [181 ] . in neutron rich matter, relevant for compact 
star cores, the theoretical predictions differ amongst themselves by more than 100 per cent 
for basic quantities such as the pressure at n ~ n sat 

53- 

Extrapolating to higher densities 



further increases these differences, on top of which new phenomena such as pion and kaon 



condensation 20|, |2l| and the presence of hyperons (such as A, £ ) only add to enlarge the 



uncertainties in the EoS. 

At some critical density, nuclear matter is expected to undergo a phase transition to 
deconfined quark matter, which is theoretically well understood only at asymptotically high 
densities, where the QCD coupling a s is small j22|. There, the stable ground state of quark 
matter is that of a color superconductor, but it is not known whether such a state persists to 
densities closer to the deconfinement transition, or whether normal unpaired quark matter, 
or some novel phase, becomes favored. Even the possibility of the normal quark phase being 
the fundamental ground state of nuclear matter (with atomic nuclei b eing only metastable) 



has been suggested in the so-called strange quark matter hypothesis [23H25| . This opened 
up the possibility for entire stars being made up of self-bound quark matter ('strange stars') 



Interestingly, despite all the advances in our understanding of QCD, most of the analysis 
of cold but dense quark matter still continues to be performed using the MIT bag model 
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dating back 35 years 27|. In this model, the interactions between quarks are absorbed into 
a phenomenological 'bag constant', which is not calculable within the model but effectively 
generated by the QCD interactions (see Ref. [28j for a discussion of this issue), and is simply 
added to the pressure of a non-interacting system. We believe that a refinement of this 
model, using a perturbative EoS for quark matter evaluated with a running a s and strange 
quark mass, should be quite superior to the plain bag model, and — absent advances in 
truly non-perturbative methods — should replace the latter whenever aiming for at least 
semi- quantitative results. 

To this end, in this paper we consider the perturbative evaluation of the QCD pressure at 
zero temperature, where the state-of-the-art result is still the pioneering order a 2 calculation 
of Freedman and McLerran (29I . 30] and Baluni j3l[. These authors, however, only included 
effects of the strange quark mass up to order a s , dropping the mass entirely at order a 2 . As 
present day knowledge suggests a strange quark mass of about 100 MeV 1321, with atomic 
nuclei corresponding to a quark chemical potential of roughly 300 MeV, one can expect 
non-negligible strange quark mass effects in the EoS (cf. Ref. |33|). In view of the situation, 
we believe that a perturbative calculation of the cold QCD EoS to order a 2 — including the 
complete strange quark mass effects — is long overdue. This provides the motivation for us 
to take on this challenge in the present work. 

Our paper is organized as follows. In Section II, we introduce our notation, explain 
how renormalization is performed, and outline the general structure of the computation. 
In Section III, we then go through all the different parts of the calculation, presenting 
the results for the individual terms and in the end assembling the entire grand canonical 
potential of the system. Section IV is devoted to a detailed analysis of our result, covering 
aspects such as the choice of the renormalization scale and the dependence of the result on 
the strange quark mass. Having gained control of the perturbative EoS, in Section V we 
consider various applications of it, studying the scenarios of stable strange quark matter and 
a phase transition between ordinary quark matter and the hadronic phase. In Section VI, 
we finally consider the implications of our work on astrophysical systems, while in Section 
VII we draw our conclusions. Several technical details, as well as most of the partial results 
of our computation, are left to Appendices A-E. 



II. SETUP 



The equation of state of a thermodynamic system is dictated by the functional relation 
between some fundamental quantity, such as the pressure or energy density, and various 
(usually intensive) parameters, such as the temperature and different chemical potentials. 
In the grand canonical ensemble, it can be solved from the grand potential, or Landau free 
energy, 

= E-fiN = -Tin Z = -PV, (1) 

where E is the (microcanonical) energy, and Z the partition function of the system. In this 
paper, we set out to perform a perturbative evaluation of the grand potential of QCD to order 
g 4 = (47ra s ) 2 in the strong coupling constant, keeping the temperature at zero but assuming 
the quark chemical potentials to be high enough so that, due to asymptotic freedom, the 
expansion converges to a satisfactory degree. Various thermodynamic quantities can then 
be obtained from the grand potential, which itself is determined by computing (minus) the 
sum of all the one-particle-irreducible (1PI) vacuum graphs of the theory. 
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Due to the infrared (IR) sensitivity of Q, at order g A it is no longer sufficient to only 
consider the strict loop expansion of the grand potential. In addition, the so-called plasmon 
sum must be performed to its full leading order, which implies resumming all gluonic ring 
diagrams containing an arbitrary number of insertions of the one-loop gluon polarization 
tensor. Upon performing renormalization, the expansion of the grand potential can be 
written as the IR and ultraviolet (UV) finite sum of (at least semi-) analytically obtainable 
one-, two- and three-loop two-gluon-irreducible (2GI) diagrams, as well as a numerically 
computable plasmon integral. 

The aim of this section is to introduce our notation, explain how the theory is renor- 
malized, and demonstrate how the computation can be divided into several distinct pieces, 
which will be addressed one by one in Section III. 

A. Notation and Conventions 

Let us consider a system of N t flavors of massless quarks and gluons, to which we add 
one massive flavor with a renormalized mass m, having in mind Ni = 2 and m = m s , the 
mass of the strange quark. We use the Euclidean metric, g^ v = 5^, and consequently define 
QCD (in the presence of finite quark number chemical potentials through the action 
Sqcd = Jd d xC QC D, with 

£qcd = \f; u F; u + M^D» + - ^ 7o )^. (2) 

As usual, the subscript B denotes bare quantities, and we have defined 

F; u = d,A a v - d v Al + g B f abc AlA c u , 

Ap = A a ^T a . 

In accordance with the metric used, our gamma matrices are all Hermitian, 7 M = 7^, and 
satisfy {7^,7^} = 2<5 M „. The flavor index % runs from 1 to Nf = Ni + 1, with = for 
1 < i < Ni and m,N f — m. Throughout the paper, dimensional regularization in d = 4 — 2e 
dimensions will be used to regulate divergent integrals, and hence /i — 0, . . . , d — 1. As 
working at zero temperature implies dealing with d = 4 — 2e dimensional integrals rather 
than discrete sum-integrals, the integration measure becomes 




where 7^ = 0.5772 ... is the Euler-Mascheroni constant, and the renormalization scales A 
and A of the MS and MS schemes have been introduced for later convenience. The chemical 
potential is accommodated in the Feynman rules by shifting the zeroth components of the 
fermionic momenta by p — >■ p + ifa, where fa is the chemical potential of the quark flavor 
in question. We choose to denote the chemical potential corresponding to the massive quark 
by HN f = fi- 

Our notation for the Feynman diagrams is as follows: Solid lines with an arrow correspond 
to quarks and wiggly lines to gluons. We work in the Feynman gauge throughout the 
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calculation, and hence our propagators have the form 

1 —if + m 

if + m P 2 + m 2 



wwvw 



P 2 



consistent with the choice to employ the Euclidean metric. Here and in the following, capital 
letters are used denote Euclidean four- momenta such as P = (po, p), with P 2 = pi + p 2 . We 
will also frequently use the abbreviations 



E(p) = Jm 2 + p 2 J u = Jj2 2 — m 2 } z = u — rh In 



l+u 



771 



(3) 



where we have introduced the dimensionless parameters rh = m/ fi,u = u/ [l. In addition, 
the group theory factors that appear in the calculation read 

d A = 5 aa = N 2 - 1, 
C A 5 cd = f abc f abd = N c 5 cd , 

N 2 -l 

C F 5ij = {T a T a )ij = — ^ — 8ij, 

where the number of colors iV c will later be set to three, and we have normalized the 
fundamental representation generators according to Tr T a T h = \5 ab . 

Finally, all expressions in this work are given using natural units, in which h = c = ks = 1, 
except when explicitly stated otherwise. 



B. Renormalization 

We perform all our calculations in terms of the bare parameters appearing in the La- 
grangian of Eq. fl2]), and only in the end express them in terms of the physical, renormalized 
ones. By doing this, one avoids dealing with explicit counter terms, but on the other 
hand cannot set e to zero until the end of the computation. The relations between the 
bare and renormalized mass and coupling constant can be obtained from the literature (see 
e.g. Ref. 0), 

m B = Z m m 
(9b) 2 = Z g g 2 

where we have defined 

Si = -Toe" 1 6 2 = C F (fC A + \C F - N f ) e- 2 - fe' 1 , 6 3 = -f3 e'\ 

with 0o and ji being constants that will be given below. As noted above, all of our calcu- 
lations are performed in the Feynman gauge, which is defined by the unrenormalized gauge 
parameter £g taking the value £b = 1- The gauge parameter renormalization constant then 



1 + 5! 



(Air) 



+ 5 2 



(An) 



+ 0(g 6 ) m, 



l + 5 3 



(Anf 



+ O(g') g* 



(4) 
(5) 
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becomes equal to that of the gauge field wave function, but neither of them enters the calcu- 
lation. The same applies to the fermion wave function renormalization constant, and hence 
we refrain from quoting the corresponding results here. 

Closely related to the above quantities are the renormalization group equations for the 
finite renormalized parameters g and m, which — given to higher order than would be 
required by our computation — read [32 



dg 2 



d In A 2 

dm 
d In A 2 



-pc 



(4 



7T 



2ft 



gr 



(4tt) 



4/3 2 



(4tt) 



-m 



7o 



(47T) 



+ 7i 



(4tt) 



+ 72 



(4tt) 



(6) 
(7) 



All constants appearing here are available from Refs. |35H38| and have the forms 



PO = 3 , Pi 

^ 2 = W + I ^f^/ ~ W C A C F N f - C\N f + || C F N] + ^ CyVj, 



7o 



3C f 



7i 



C f 



07 



-AT f 

3 / 



72 



CV (±f C 2 - if + ij^Ci + C F N F (-23 + 24C(3)) 



/•' 4 



-f - 24C(3)) - fiV 2 



(8) 



The renormalization group equations may furthermore be integrated to give the running 
coupling constant a s (A) = g 2 (A)/(47r) and strange quark mass m(A) as functions of the 
renormalization scale. Following the treatment of Ref. [38j by expanding the results in 
powers of a s , we obtain 



a s {A) 
m(A) 




L = In (AVA^ 



70 //9o 



1 + Ai MA) + A^M 



(2GeV) / 



1 + A x 



q s (2GeV) . A|+A 2 ( a s (2 GeV) 



+ 



^ a s (2GeV) ^ 



2 ' 



(9) 
(10) 



with 



A x 
A 2 



Prya _7i_ 

2/3 2 + 4/V 




/3i7i 



+ 



72 



A 2 160o 



(11) 



Here, A^g- denotes the MS renormalization point, and we have chosen a fiducial scale of 2 
GeV where the strange quark mass takes the value m(2 GeV) ~ 100 ± 30 MeV j32| . 

One way to fix the renormalization point is by requiring that the coupling of Eq. ([9]) agrees 
with the current best experimental value for the parameter at the fiducial scale, ct s (2 GeV) = 
0.2994 + 0.0152 - 0.0141 , leading to A^ = 0.378 + 0.034 - 0.032 GeV. For completeness, 
we note that an alternative to this would be to use a lattice QCD determination of the 
deconfinement temperature at small baryonic densities. In one lattice study, this was found 
to give Ams ~ T c /0.49 [39], which — depending on the result for T c one uses [40J, |41 



8 



(a) 




I , I , I , I , I , I , I , I , I , l_ 

0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 

A [GeV] 



(b) 




0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 
A [GeV] 



FIG. 1: Renormalization scale dependence of a s and m from Eqs. ([9]), (|10p . using values of A 
relevant for compact stars. The full lines and dark shaded regions correspond to the central value 
and uncertainty of A^g = 0.378 + 0.034 — 0.032 GeV, respectively. The light shaded region in m 
corresponds to the uncertainty in m(2 GeV). The circles with error bars correspond to reference 
points obtained from Ref. 32J. 



results in A^g ~ 0.375 ±0.016 GeV. As this result happens to lie within the uncertainties of 
the first method, we thus adopt A^g" = 0.378 + 0.034 — 0.032 GeV for the remainder of this 
work. 

In Fig. [U we display the functions a s (A) and m(A), as given by Eqs. © and f|T0|) . to 
provide an indication of the value and uncertainties of the quantities. From here, it is 
evident that below A ~ 0.8 GeV, the results become dominated by the various uncertainties 
discussed above, setting a lower limit for the values of A for which one can even optimistically 
expect quantitative results from any perturbative calculation (even in the hypothetical case 
that one was able to determine the quantity in question to infinitely many orders in o^). 1 
Reducing these uncertainties would require a better determination of both a s and m at the 
fiducial scale, which is not the aim of our work. 



C. Organizing the Calculation 

To the perturbative order we are working in, the grand potential of QCD obtains contri- 
butions from the (renormalized) single quark loop fin,, the two- and three-loop two-gluon 
irreducible (2GI) vacuum diagrams of the theory f2 2 ci, as well as the plasmon ring sum n p i as , 

^ = ^1L + ^2GI + ^plas- (12) 



1 In principle, one should adjust Nf in Eq. (jHJ when crossing the strange quark mass threshold (cf. Ref. [42 1 
for a pedagogical review on heavy quark decoupling). However, from Eq. (|10|) one finds that the decoupling 
of the strange quark occurs at A ~ 0.6 GeV, where our results are not quantitatively correct in any case. 
Hence, we will simply ignore this effect in the following. In contrast, note that we will find that setting 
r7i = would change our results for the partition function at A > 0.6 GeV, so the strange quark mass 
effect cannot be altogether ignored. 






.0 .,.,.,.„. 

FIG. 2: The two- and three-loop 2GI diagrams contributing to the grand potential of QCD. 

The 2GI graphs are displayed in Fig. [2J while the plasmon sum f2 p i as corresponds to the 
sum of all gluonic ring diagrams containing at least two insertions of the one-loop gluon 
polarization tensor, depicted in Fig. [3^ (c/. Appendix [X]) . As explained in more detail in 
Section UlIDI as well as Ref. jsj], the benign IR behavior of the vacuum (T = fi = 0) part of 
the tensor allows the division of the ring sum into three parts. These include the three-loop 
vacuum-vacuum (VV) and vacuum-matter (VM) graphs of Fig. [3}o and c, as well as the 
'matter' (vacuum subtracted) ring sum of Fig. [3jT, 

^plas = ^VV + ^VM + ^ring- (13) 

Here, the VV diagram only produces a chemical potential independent contribution to the 
grand potential, and is thus neglected. The matter part of the polarization tensor merely 
includes the fermion loop, which behaves like 1/P 2 in the UV, so one finds that the ring 
sum Q T i ng is both IR and UV finite (c/. Section UlIDp . 

To further organize the calculation, one may separate the contributions of the massless 
quarks to both f2 1L and f2 2 Gi from the rest by writing 

n 1L = n r {i=° + n?L, ^iT = n 1L (n = m), (14) 

J1 2 GI = ^2GI° + ^2GI> ^2GI° = ^2Gl(.H = m ) , (15) 

and defining fi^L an d ^2Gi as the respective differences. With the vacuum-matter diagram, 
the issue is slightly more subtle, and we in fact have to write it in three parts, 

(16) 

Here, ^vm * s defined as the VM diagram composed of the massless part of the matter 
polarization tensor and the m — > limit of the vacuum polarization tensor (i.e. containing 
contributions from Nf = Ni + 1 massless quarks), while fiyM consists of a massive matter loop 
and the entire, m dependent vacuum polarization tensor. Finally, Oy M is the contribution to 
the grand potential originating from a diagram composed of the massless part of the matter 
polarization tensor coupled to the difference of a massive and massless vacuum quark loop. 
It is easy to verify that the sum of these three functions indeed equals Qvm- 

Collecting the various pieces of the above expansion, the result for the grand potential 
may be written in the form 



n = n m =° + n m + n x VM + n ring , (17) 



where we have denoted 



rr=° = fir=° + fi™=° + ny=°, (is) 

fl m = + ^2GI + ^VM- (19) 
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FIG. 3: The generic form of the ring or plasmon sum. 



The reason for the above definitions now becomes clear: Having written f2 in terms of alto- 
gether eight different parts, we have managed to separate the contributions of the massless 
quark flavors (up to their contribution to the ring sum) into the function Q m=0 , the value of 
which may be directly taken from Ref. In the massive part, the contributions from the 
single quark loop Q 7 ^ (including renormalization corrections to it), the sum of the two- and 
three- loop 2GI diagrams ^2Gt, as wen as the VM part 0™ M have also been separated. In the 
following section, we will proceed to evaluate these various parts of the grand potential one 
by one. 



III. THE COMPUTATION 

This section contains a discussion of the evaluation of each of the independent contribu- 
tions to the grand potential, as defined in the 'master formula' of Eq. f TTTj) . For brevity of 
presentation, many of the calculational details as well as nearly all intermediate results are 
left to be listed in the Appendices. Readers not interested in the calculational details may 
wish to skip this section altogether and proceed straight to Section HV] where the final result 
is analyzed. 



A. n m -°: Massless Quarks 



The contribution of the massless quark flavors to the grand potential, including the 
massless single quark loops fi^ =0 , the sum of the two- and three- loop 2GI graphs containing 
massless quarks fi^ , and the massless VM graph fiyj^ can be extracted directly from 
Ref. §. One finds 



o™.=o 

"lL 

V 

Ctm=0 
"2GI 

V 



Oin=0 
"VM 



V 



127T 2 f/*' 



4tt 2 
2 



47T 



i=l 



5C A - 2N f 8 



3e 



Ca + H Cf -™n 

3 2 3 ' 



ff 2 (A) 
(4vr) 2 



N, 



ATI 2 



((5/2 + e)C A -N f )J2rf 



i=l 



2 , A 
h 4 In 

3e 2 ft 



Nf) In 



A 

2/ii 



52 \ g\A) 
9 J (4tt) 4 



(20) 



(21) 
(22) 
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where V is the volume of the system, renormalization corrections have been taken into 
account, and the result is expressed in terms of the renormalized gauge coupling constant g. 
Adding the three pieces together, one obtains for the function Jl m=0 



n 



m=0 



N, 



V 



+ d A 



4tt 2 
4 



4 J N c g\K) 



Nf 



11C, 



In 



A 

2 fa 



142 ^ r7 22 Ar 1 - 4 

9 2 9 ' 



g 4 (A) ) 
(4vr) 4 J' 



(23) 



B. Q m : Massive Quark 

In this subsection, the contribution of the massive quark to the grand potential of QCD 
is evaluated by separately considering all the parts of the function Q m , as listed in Eq. ( fl9l) . 
The final result for the function will be given at the end. 



a. The Massive Quark Loop 

The grand potential of one free massive quark flavor at zero temperature but finite chemical 
potential can be written in the form 



n 



1L 



V 



-2N r A 2e 



d 3 " 



2c 



V 



{h-E{ P ))B{h-E{ P )), 



(24) 



(27r) 3 " 2e 

where the regulator e has been kept non-zero because of the bare mass parameter appearing 
inside of E(p). Writing the integral in the form 



n 



1L 



V 



c (2vr) 3 - 2e 



\l up— mi. 

V B d PP 2 - 2 * 



V ~ VP + W| , 



(25) 



with S n denoting the area of an n-sphere, one can perform a power series expansion in g 
and e. Keeping terms up to and including 0(g 4 ) and 0(e 2 ), the integral is evaluated to be 

d 3 p 



"lL 

V 

— Sim 



-2N r A 2e < 



(27T) 



3 (// - E(p)) 9(h - E(p)) 



(26) 



2 ^ 2 (A) 



(4tt) 



5\m 2 



S 2 m 2 + 



l + 2m 2 







dm 2 



/(A) \ f d 3 - 2 -p OJji - E(p)) \ 
(4vr) 4 J J (2vr) 3 - 2e E(p) J' 



in which the mass and coupling constant are now the renormalized ones. 
A straightforward evaluation of the remaining integrals finally gives 



n 



1L 



V 



12tt 2 

,2 



j • • > . 3/// 

uji\ji — -m JH ; — hi 



H + u 



(4tt) 2 



+ d A m < C A 



11 
72" 



- AT 



"2 5 " 











97 
6e 

ff 4 (A) 
(4tt) 4 



/i + C f 



18 3 

— - - \I X + 18/ift 



(27) 



where the I n (m)'s denote integral functions defined in Appendix Here and in the follow- 
ing, we have suppressed the arguments of these functions, except when they are necessary 
to avoid confusion. 
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b. The Massive Two- and Three-loop Skeleton Graphs 
The four massive 2GI graphs of Fig. [2] contribute to the grand potential as 

- = D a + D b + D c + D d , (28) 

where the functions Di appearing on the right hand side correspond to the diagrams a-d. 
Their evaluation proceeds in several steps. First, we write down expressions for the graphs in 
terms of the Feynman rules and contract all color and Lorentz indices, reducing the original 
functions to sums of scalar integrals. The fmiteness of the chemical potential prevents the 
application of integration-by-parts identities at this point, so instead we perform the p 
energy integrations using the residue theorem. Collecting the outcome of this step, one 
observes that each graph hereby reduces to a simple form that can be viewed as the result 
of having performed all possible 'cuts' on the fermionic lines of the diagram, their number 
ranging from zero to the number of loops. 

Each time a line is cut, the corresponding propagator is placed on shell and the remaining 
part of the diagram integrated with respect to the three-momentum in question with the 
weight 

d 3 p 9{iM-E(jp)) 



(2tt) 3 2E(p) ' 

while /i is set to zero in all other propagators. This implies that the integrands of the phase 
space integrals become d = 4 — 2e dimensional vacuum (T = /x = 0) two-, four- and six-point 
functions, which may be evaluated using integration-by-parts relations, conveniently imple- 



mented e.g. in the Mathematica package FIRE 43] and results available in the literature 
[see e.g. Refs. HJ447I)- 



In the end, the cutting procedure separates each n-loop graph into n + 1 parts 

A = A° c + D} c + Df + Df + • • • , (30) 

corresponding to the number of cuts performed. Thus the entire 0(g 4 ) 2GI contribution to 
the grand potential separates into three pieces 

C\m Qm,lc i Q?n,2c . qt?i,3c /q-, \ 

"2GI — Ji 2GI ' 2GI "r"2GI- 

Here, we have taken advantage of the fact that the zero-cut contribution is always indepen- 
dent of the chemical potential and may thus be dropped. We also note that the computa- 
tionally most complicated triple-cut part of the grand potential is automatically UV finite 
and thus in no need of regularization, which simplifies the computations significantly. 

As the detailed evaluation of the function is quite lengthy and almost void of physical 
content, we leave it to Appendix [Bj There, we first go through the scalarization and cutting 
procedures for each graph, and in the end list the results for the functions Q™g\ and 

a^?j in Eqs. (JBUD (TBI9]). 

c. The VM Graph 

The evaluation of the massive vacuum-matter graph fiyM proceeds in a way highly analogous 
to that of the 2GI graphs, i.e. via first performing the Lorentz and color contractions and 
then taking one and two cuts in the remaining scalarized diagrams (the number of fermion 
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propagators limits the maximum number of cuts to two here). The result can be written in 
the form 



VM 



"VM ' 



m,2c 
VM > 



(32) 



where the functions fi™M C ) ^vm are given in Eqs. (1C4j) and ( 1C5|) . 



d. Assembling the Result for fl m 

Collecting the results of this section and Appendices |B] and [Cj we see that the contribution 
of the massive quark flavor to the grand potential is of the form 



Mi + M; 



ff 2 (A) 
(An) 2 



M. 



g 4 (A) 
(4tt) 4 



(33) 



where the two first terms of the result have been calculated before, e.g. in Ref. 33]. In 



our presentation, the function Mi, corresponding to the grand potential of one free massive 
quark flavor, can be read off from Eq. (J2j 



Mi 



NcfS 
2471- 2 



2u A - 3zm< 



(34) 



where we recall the abbreviations introduced in Eq. ([3]). At the next order 0(g 2 ), the 
function Q m on the other hand obtains contributions from the two-loop graph of Fig. 
and the 0(g 2 ) renormalization correction to the massive single quark loop. Adding these 
functions together using Eqs. ( 1271) . ( IB 171) and ( 1B18j) . we find 



M 2 = -6zm 2 ln- + 2u 4 -4zm 2 -3z 2 

An 2 m 



(35) 



To obtain the 0(g 4 ) contribution to Q m , one needs to add together the corresponding 
parts of the functions fi^-j and Oy M . Due to the length of the resulting expressions, 
we organize the calculation in terms of the number of cuts and write Ms in the form 



M, 



M\ c + M\ c + M 



3 • 



(36) 



dealing with the three pieces separately. 

There are three different types of single cut contributions to M^. The order g 4 renor- 
malization correction to Q 7 ^ from the last two lines of Eq. (I271) . the function VL™qi from 
Eq. (1B17I) . as well as the single cut contribution to the massive VM graph f2y M from Eq. (1C4I) . 
Summing up the various pieces, one finds 



M\ c 



(2n) 2 
3 



- C F 
+ N, 
+ Am 2 C F 



A 



' . A 185 
22 In 1 

m 3 



, A 1111 
In 1 

m 24 



4vr z 



+ 4vr 2 ln2-6C(3) 



12 In — + 5 In — + 



A 313 35tt 2 



+ 



m 



A 



m 



A 71 



- 6 In — + 13 In — + — + 



6 

2tt 2 



in 



12 



-87r 2 ln 2 + 12C(3) 
+ 6-2tt 2 )m 2 z 



A 



A 



3 3 In — +4 In- + 4 



rn 



rn 



u 



(37) 
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where it is interesting to note that all the 1/e divergences have canceled. This implies that 
no cancellations need to take place between the single and double cut contributions, which 
is an important computational simplification. 

At the two-cut level, there are two sources of contributions to The two- and three- 
loop 2GI diagrams from Eq. (IB 18 1) and the two-cut piece of the function fly M from Eq. ( 1C5I) . 
The sum of these parts reads 



M 



2c 



16 



62 



10 



22 



j—^ = d A \c A y - —it + —m fh + g/ic - Y m ^ + J io - y [A - 2m J 2. 



m , 



+ C F {In + [2A(m 2 I 2 - m 2 I lb h + 2m A I 2b ) + 48m 4 J 8 

"4 



In 



A' 



m 



Ar /10 t2 20 2 T 2 T 4 2t 
+ N f —I 2 - —m 2 I 2 - -I lc + -m 2 I 2c + 



9 



9 



-I 2 - -m 2 I 2 



In A )-§ 

ml 3 



(38) 



where the definitions and values of the various integrals I n are given in Appendix [D] 

Finally, the only triple cut contribution to the grand potential originates from the 2GI 
graphs of Figs. [2b, |2b and[2ji, implying that it can be read off directly from Eq. (IB19j) . One 
finds 



Mf 



-d A {C A 



2IJ 2 - 4J 5 + 8m 4 / 6 - 4/ 7 



+ C t 



2l\l xh - AI X I 2 - 8m 2 hl 2b + 8m 2 I 3 + 8m 4 / ; 



3b 



2/4 + 8/5 - 16m 4 J 6 + 8J 7 - 8m 2 /i/ 8 + 8m 4 / 9 



(39) 



Summing up the single, double and triple cut contributions to A4 3 leads to a complicated 
expression, in which several parts must be evaluated numerically. To this end, we have 
found it most convenient to express the result in terms of a basis of simple fitting functions, 
leading to the result 



M, 



d A fi 4 



2tt 2 



m 



[UC A - 2N f ) z + 18CV (2z - u) 



In- 

m 



+ 



3 

'/ 

+ C A 



C A 22u 4 



185 



zm 



9C 



33z 2 J + — - (l6m 2 u(l -it)- 3(7m 2 - 8u)z - 24z' 



- N f (Au A - 13,2 m 2 - 6z 2 



in- 

m 



11 , m 71 „ , A A „ /17 „ 
y In y - y + 0i(m) J + C F ( — + 2 {rh) 



+ JV> ( - In — + — + g 3 {m) J + g A {m) \, 



(40) 



,3 2 9 

where the functions Q n (m) (that are defined so that they vanish in the m — > limit) read 

& (m) = 327r 4 m 2 [ - 0.01863 + 0.02038 m 2 - 0.03900 m 2 log(m) 

+ 0.02581 m 2 (log(m)) 2 - 0.03153 m 2 (log(m)) 3 + 0.01151 m 2 (log(m)) 4 1 , (41) 
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02 (m) 



03 (m) 



04 (m) 



= 327r 4 m 2 - 0.1998 - 0.04797 log(m) + 0.1988 m 2 - 0.3569 m 2 log(m 
+ 0.3043 m 2 (log(m)) 2 - 0.1611 m 2 (log(m)) 3 + 0.09791 m 2 (log(m)) 4 " 
= 327r 4 m 2 [ - 0.05741 - 0.02679 log(m) - 0.002828 (log(m)) 2 
+ 0.05716 m 2 - 0.08777 m 2 log(m) + 0.0666 m 2 (log(m)) 2 

- 0.02381 m 2 (log(m)) 3 + 0.01384 m 2 (log(m)) 4 ] , 

= 32vr 4 m 2 [0.07823 + 0.0388 log(m) + 0.004873 (log(m)) 2 

- 0.07822 m 2 + 0.1183 m 2 log(m) - 0.08755 m 2 (log(m)) 2 
+ 0.03293m 2 (log(m)) 3 - 0.01644m 2 (log(m)) 4 ~ 



(42) 



(43) 



(44) 



C. ^vm : The Cross Term VM graph 

Next, we look at the vacuum-matter term f2y M , in which the quarks in the matter part of 
the diagram are massless, while the vacuum part corresponds to the difference of a massive 
and a massless quark loop. Performing again one and two cuts and using the results of 
Appendix [UJ one obtains 



QX 

V 



di 



g 



m 



(4vr) 2 12tt 4 



E 



x 



dw In 



m 2 w 
Pi 



w 



w 



-arctanhU- 



w 



w 



+ w log(4u>) 



d, 



g 



m 



4 N t 



5> 



fJLi 



(An) 2 12vr 4 ^ \m + ^)' 
where the integral can be approximated by the pocket formula 



(45) 



Ut) 



-3t 4 (1 - In t) 



0.83 



+ 



0.06 



0.056 



{l-tf (1-t) 

+ y ~ ll f 1 - 005 - 0.272*(1 - t) + 0.154t(l - tf 
t{ 1 — tr v 



(46) 



D. ^Hng : The Plasmon Contribution 



The setup for the plasmon sum calculation can be found in standard textbooks, e.g. Sec- 
tion V of Ref. [11], according to which the contribution can be written in the form 

^plas 



V 

(U f d A K 

T 



(2vr) 4 L 
tes t] 

G(k ,k) 



2 In 1 



G(k Q , k)\ + ]n ( 1 _ F(h,k) \ + 2G(k ,k) + F(ko, k) 



K 2 J \ K 2 J K 2 K 2 

where k denotes the magnitude of the three- vector k. In this equation, we have defined 

1 



(47) 



2 n»(k ,k) - ^u 00 (k ,k)\ , F(k ,k) = ^n 00 ^,*;), 



(48) 
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with 11^ = IIq — II-, where IP^ is the one-loop gluon polarization tensor (cf. Appendix [A]) . 
One should note that we have used the notation F, G, rather than A 1; A 2 sometimes found 
in the literature 30|, |49) to avoid confusion with the renormalization scale A. 

Taking advantage of the benign IR behavior of the vacuum polarization tensor, we ex- 
pand the logarithms in powers of the vacuum tensor, keeping only contributions up to and 
including 0(g 4 ). This leads to 



In 1 



In 1 



G(k ,k) \ 
K 2 ) 

G ma t(^o ? k) 
K 2 



G vac (k , k) G vac (fc , k)G mil t(k , k) 1 G 2 WSJC (k , k) 



K 2 



K 4 



K A 



(49) 



and a similar expression for the F(k , k) part. This separation allows us to split the plasmon 
sum into the VV, VM and ring sum pieces, the first two of which were already covered in 
the previous section. The last part, which is considered here, reads 



firing 



2 

+ In (1 



d A K 



2 In - 

(k ,k)\ 



K 2 



<Jmat(^0) k) \ 

K 2 J 

2G nmt (kp, k) 
K 2 



+ 



F mat (k ,k) 
K 2 



(50) 



To isolate the g A In g term from the rest of the plasmon sum, we separate from Eq. ([50 
the combination (both for F and G) 



d A K 
(2tt)* 



In 1 - 



F mat (K = 0,$Y 
K 2 



+ 



F mat (^ = 0,$) , F2 at (^ = 0,$) 



K 2 



+ 



2K 2 (K 2 + X 2 ) 



(51) 



where \ is a fictitious mass scale that will drop out in the final result and we have switched 

k 



from using the variables (fco, k) to K = y k$ + k 2 and $ = arctan^. The reason for 
introducing this particular function is that it captures the relevant IR physics in a form 
where the K integration may be performed analytically, giving 2 



(2*) 



rf$ sin 2 $ — 



-i + ln(-F mat (ir = 0,$)/r 



(52) 



The complete ring sum contribution to the grand potential then becomes 



firing 



ring 1 J "rmg 1 "-""ring; 



)(2) 



)(3) 



(53) 



where the 0(g A In g) contribution has been isolated in 



fi 



(i) 

ring 



V 



d A 9 4 \ng ^l 2 
(2tt) 3 



rf$ sin $ 



2 Gg Mt (X = 0,$) | FlJK 



0,$) 



r 



(54) 



2 Note that there is an error in the corresponding Eq. (5.5.2) of Rcf. 



48] 
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and the remaining part of the g 4 contribution is contained in the functions 



n 



(2) 
ring 



V 



d-A g A 
' 2 64tt 7 Jo 



tt/2 



d$ sin $ 



-2n 2 F mat (K = 0Y 



1 + ln 



4 4 



F mat (^ = o) 



1 / -47T 2 G mat (ir = 0) - 

2 I 



1 + ln 



4 4 

XT 



GLt(^ = 0), 



(3) 
ring 



V 



d 4^rdKKr /2 d* sin 2 $l 
2 (27r) 3 io Jo g 4 



Fl at (K = 0)-F 1 



mat 



K 2 



FLt(K = 0) + 2^ at (K = 0) - 2GLt 2GL t (^ = 0) 



K 2 + X 2 



K 2 



K 2 + X 2 



(55) 



Collecting all of the above formulae, one observes that the result for the plasmon ring 
sum can be written in the form 



^ring 



2 In 



d A g A f ( 

5WlH 

1 / 19 2tt 2 J 15 (/2) 16 



9_ 

Arc 



2 [ 3 4 3 



f Tjffi + y - ln2 ) In 2 + J 16 ( m, /} 



1=1 



( 14 



2 In 



47T 



1 / 16 / -2 

- fl 17 (m,/l«) + y(l - ln2)ln 2 J 18 + 1 19 fm,/i 



(56) 



+ f 



'13 



2 In 



47T 



1 /' X6 

2 ( J 20 + y (1 - ln 2) ln 2 I 21 + J 22 ( m, /T 



where /2 = // 2 , • • • , /MrJ and // = while the definitions of the J n 's — including 

approximation formulas to the numerically evaluated functions — can again be found from 
Appendix [Ql Finally, we note that upon sending // — > m in Eq. ( 156|) and using the result that 

in this limit I\q fm, /t ^ — )■ —0.85638 . . ., the result of iVj massless quark flavors presented 

in Ref. jEf is recovered. Similarly, one can show that taking m — > 0, one obtains the correct 
result for iVj = Ni + 1 massless flavors. 



E. The Result for the QCD Grand Potential 



According to our 'master equation' ( 1T7|) . the QCD grand potential, evaluated to the 
perturbative order g 4 in the MS scheme, can be written as the sum of four terms, each of 
which we have evaluated in the previous sections. Of these functions, 

• Q m=0 corresponds to the contribution of the massless quark flavors, and is available 
from Eq. (|23l). 

• Q m corresponds to the contribution of the massive quark, and is available from Eq. ( |33jl . 
using Eqs. (1341), ([35]) and (j40|). 
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• ^vm corresponds to a cross term vacuum-matter diagram, and is available from 
Eq. (|25}. 

• ^ring corresponds to the plasmon ring sum, and is available from Eq. ( )56l) . using 
integrals listed in Appendix ID] 

Before proceeding to use our result to study various physics questions, we want to stress 
that it has passed all the consistency checks available in the literature. In particular, the 
result is renormalization group invariant, i.e. the explicit logarithms of A match with those 
originating from the running of the coupling and the strange quark mass, and the result 
correctly approaches the cases of two and three massless quark flavors in the respective 
limits ji — > m and /i/m — > oo. Of these, the fact that the massless Nf = 3 result of 
Refs. [29| and j^J is analytically reproduced is highly non-trivial, and a strong indication 
that our calculation has been performed correctly. 

At this point, a remark regarding the region of applicability of our result is in order. As 
can be verified from Eqs. f)37p and (140 p . the grand potential contains a part that behaves 
in the ft — )■ m limit as g 4 fj 2 m y /fi — m, leading to a divergent 1 / -y/// — m contribution to 
the strange quark number density. Inspecting the lower order contributions to Q, one can 
identify the reason for this behavior: In one part of the result, the expansion parameter is 
g 2 fi/(fi — m), which diverges in the limit m — > fi. It would be interesting to analyze this 
effect in more detail and investigate what type of a resummation scheme would be required 
to describe the limit properly, but this is beyond the scope of the present work. For the 
purposes of this article, we merely note that as long as one confines the analysis to values 
of \i satisfying 

> g\ (57) 

m 

no problems will occur. Due to the rather low value of the strange quark mass, this in 
practice provides no extra limitation for the applicability of our result. 



IV. ANALYZING THE RESULT: QUARK MATTER EOS 

In this section, the properties of the cold quark matter EoS are derived from the result 
presented in the previous section, and issues such as the strange quark mass dependence 
and the choice of renormalization scale will be discussed. We want to stress that while 
the calculation of the QCD grand potential is completely unambiguous, the extraction of 
physical quantities such as the EoS makes it necessary to adopt particular strategies for the 
use of the result, e.g. regarding to which quantities to truncate at a given perturbative order 
while still ensuring thermodynamic consistency. With this caveat in mind, it is pleasing that 
many properties of the resulting EoS are naturally meaningful and intuitive. For example, 
in this section it will be demonstrated that the result for total quark number density for two 
massless and one massive quark flavors smoothly interpolates between the cases of two and 
three massless quarks. 



A. Choice of Renormalization Scale 



With the functional form of a s (A) and the strange quark mass m(A) specified by Eqs. (Q 
and (fTOl) . the only undetermined parameter in the perturbative EoS is the choice of the 
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renormalization scale A. At high temperature and low density, the canonical choice for it is 
the first non-vanishing Matsubara frequency 2ttT, with which resummed perturbation theory 
has be seen to provide a reasonable description of the lattice results for the EoS and entropy 
density for temperatures at least a few times the critical temperature of the deconfinement 
transition jH, [Ho|. For one massless quark with chemical potential fi, phenomenological 



models suggest the choice A = 2ir^T 2 + (fi/ir) 2 up to an overall factor of two |5lH54|. which 
also happens to be the scale appearing in the leading order fermionic quasiparticle mass 50 . 
Finally, a comparison of the perturbative and exact pressure at zero temperature, known in 
the limit of a large flavor number, also suggests A = 2\i [55|. 

For the above reasons, our canonical choice for the renormalization scale will be A = 
2Eft/%i around which we will vary the parameter by a factor of two. We have also 
investigated more complicated choices of A taking into account the finite value of m s , but 
have found this to have only minor effects on the results. 



B. Warmup: Massless Quarks 



To illustrate our approach in a simple setting, let us first consider the case of three 
massless quarks (Nf — iVj — 3) at the same chemical potential (fi u — ^ — ^t s = fi), 
corresponding to a system that is locally electrically neutral. Up to 0(a 2 ), the perturbative 
grand potential is then given by the sum of Eqs. (1231) and (1561) . with the massive quark's 
chemical potential set to zero in the plasmon integrals. In principle, one could extract the 
pressure directly from these equations, via P = — Q/V. The determination of the EoS will, 
however, also require the total quark number density 



n = — — 



]_d_ 

Vd^L 



n 



m=0 



m=0 
ring 



(58) 



which, when evaluated using this formula, will receive corrections at higher order in a s for 
instance from -§^ot 2 . These higher order terms are beyond the accuracy of our calculation, 
and hence typically ill-behaved, but cannot be simply dropped because this would imply 
ndfi 7^ VdP and therefore ruin thermodynamic consistency. However, it is also possible 
to choose the quark number density as the fundamental quantity, keeping only terms up 
to (and including) 0(a^), and use it to determine the other parts of the EoS by requiring 
thermodynamic consistency. Using Eq. (1581 . we thus obtain (for general Nf) the three-loop 
result 



n (2) (/i,A) = n w O) 



(o), 



1-2^-f^n- 

7T V 7T / V 4 



AT , N f a s _ , A 
+ N f In + p In - 

7T u 



11 In 2 - 0.369165N/ 



(59) 



Here, (3 is given by Eq. (jSJ), and the non- interacting quark number density reads n^°\fi) = 

Nf^/Tl 2 . 

In Fig. Hk, we display the behavior of the function n^(/i, A), evaluated numerically using 
Eq. OH]) with the value of the renormalization scale varied between A = // and 4/i. To assess 
the convergence properties of the perturbative expansion, we show in the same figure also the 
two- loop result , obtained by truncating Eq. fl59|) at order a s and — to be consistent - 
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FIG. 4: The renormalization scale dependence of the massless quark number density and pressure 
for perturbation theory to order 0(a).) and 0(a 2 s ) (light and dark shaded regions, respectively), 
with Aj^g = 0.378 GeV and A ranging from \x (right boundary) to 4/i (left boundary). The dashed 
(full) lines correspond to the first (second) order results with A = 

also truncating Eq. fl9]) by setting (3i = 0. Inspecting the result, we find that the total quark 
number density n(fi, A) is (even when normalized to the free quark density) a monotonically 
increasing function of /i, and vanishes for a particular chemical potential /io(A). This critical 
chemical potential is below the regime where our computation is reliable, and the fact that 
the quark number densities become negative for fi < /i can be seen to signal the breakdown 
of our calculation. We will therefore set n(/z, A) = for fi < fio(A). 

As one can see from Fig.H^, our results for are within the uncertainty band of the two- 
loop result, given by the renormalization scale dependence of for fi > 1 GeV, indicating 
that the perturbative expansion for the quark number density converges reasonably well. 
More surprisingly, even at lower fi, where the uncertainty band generated by A//i e [1,4] 
is large, our values for /io(A) do not differ much between and n^ 2 \ This should be 
contrasted with the case of high temperatures and low densities, where the convergence 
properties of weak coupling expansions are in most cases substantially worse. 

Finally, from the quark number density one also obtains the pressure through the relation 



where B is an integration constant, equal to minus the pressure at fi = /io(A). While in 
a purely perturbative calculation B would usually be set to zero, in a realistic description 
it should be taken to be non-vanishing as the pressure can in any case only be determined 
up to an additive constant, representing the pressure difference between the physical and 
perturbative vacua. In our work, we will consider B a free parameter, which allows us to 
take into account non-perturbative effects not captured by the weak coupling expansion. In 
fact, using the free quark number density with /i = instead of in Eq. f )60|) . one 
recovers the expression for the pressure in the original MIT bag model, with B taking the 
role of the bag constant. 

When performed in the way described above, our originally purely perturbative calcula- 
tion can be seen to offer the possibility of adding non-perturbative effects to a result that 




(60) 
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is guaranteed to have the correct behavior at high energy densities. Due to physics criteria 
(e.g. requiring the energy density to be positive), the possible values for B are, however, 
typically rather restricted, allowing us to make quantitative statements that are not possible 
in the original MIT bag model. In particular, if the equation of state is to be used down to 
/i = /io, one must require B > 0, because otherwise the energy density 

e(fi,A) = -P(ji,K)+fin(p,K) (61) 

will be negative at // = fiQ. It is only if one ceases to use the perturbative EoS at some 
fi > fio — for instance when matching to a hadronic EoS at smaller fi — that B can take 
on negative values. 



C. Massive Quarks 

In a realistic attempt to describe cold deconfmed quark matter, the non-vanishing mass 
of the strange quark must be taken into account. Up to and including the perturbative 
order a 2 s1 the grand potential is then given by the main result of our calculation, collected in 
Section IIII El To extract physically meaningful information from this expression, we follow 
the same strategy as in Section IIVBI up to some minor modifications necessary to maintain 
thermodynamic consistency. 

We begin by evaluating the up (u) quark number density to 0(a 2 ), truncating terms 
of higher order, and use this to obtain expressions for the down (d) and strange (s) 
quark number densities that are thermodynamically consistent, i.e. satisfy relations such 
as dn s /d/i u = dn u /dfi s . This procedure leads to the results 

n u (Vu, Vd, AO = ^(nuj^/is), 

n d (Hu, fid, aO = / dfi' u d^ d n {2) (fi' u , fi d , (j, 8 ) + n d 2) ([i u = u , fi d , fj, 8 ), 

n s (n u , fi d , fi s ) = f d//^ s n[ 2) (/i^/i d ,// s ) + ( dfi' d d^n d 2) (u Q (^ d ^ s )^' d ^ s ) 

+nf ] (fi u = u ,n d = d , n s ) , (62) 

where the functions nf^ are defined so that they contain no terms beyond 0(a 2 ), and u 
and do are free integration functions. For reasons discussed in Section TlV Al — and to better 
facilitate comparison with the massless case — the renormalization scale is taken to be of 
the form 

Aoc^±^±^, (63) 

where the canonical choice of two for the proportionality constant will be varied by a factor 
of two to test how strongly our result depends on this choice. 

In this work, we consider deconfmed quark matter that is locally charge neutral and in 
beta equilibrium. 3 Chemical equilibrium is reached via the weak processes 

d — > u + e + v e , u + e — > d + v e , 
s — > u + e + P e , u + e — > s + v e , 
s + u <-> d + u , 



3 Relaxing the assumption of local charge neutrality to a global one in systems with a mixed phase is 
typically only a minor effect in comparison with the renormalization scale dependence [56J ■ 
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which imply the conditions 

= Vd = A* > /i u = /i - /x e , (64) 

with /z e being the electron chemical potential. The contribution of the neutrinos has been 
neglected here because for the quasi-static compact star systems we have in mind, neutrinos 
escape quickly, as their mean free path is considerably larger than the physical size of the 
system. Local charge neutrality on the other hand relies on the presence of electrons and 
leads to the relation 

2 11 

^n u - -n d - -n s - n e = , (65) 

where n e = /ig/(37r 2 ) is the electron density. Solving Eq. (I6"5"|) with Eq. (1M|) fixes the electron 
chemical potential \x e as a function of //. 

The above constraints of beta equilibrium and charge neutrality enable us to simplify the 
expressions of Eq. (1621) significantly. Namely, for a very particular choice of the integration 
functions do and uo, do(p, s ) = fi s and UoifMd, A* g ) = fi s — He(fJ>d,^s), we see that each of the 
integral terms in Eq. (162]) vanishes on the one-dimensional curve in /ij space where both 
beta equilibrium and charge neutrality are maintained. In this physical subspace — and in 
particular when calculating fM e (fj,d, /-0 from Eq. fl6"5"|) — we may use the truncated expressions 
nf 1 for all quark flavors, 

n u (/i,A) = n^\fi - /i e (/i),/x,/i) + 0(a^ 



s ) 



n d (ji,A) = n d 2 \fi - fj, e (jj),fi,,fj) + 0(a 3 s ) , 

n a (fi,A) = nf\n - /i e (/i),//,/i) + 0{otl) , (66) 

where we have emphasized the dependence of the result on A by explicitly reinstating it as 
an argument of the functions. In these results, we will again set = whenever the number 
density in question becomes negative. We also note that while in our starting point Eq. ( |62|) 
we chose n u rather than or n s as the fundamental quantity, our end result Eq. (|66|) treats 
all flavors in a symmetric fashion. 

In terms of the above results, the total quark number density is finally given by 

n(/i,A) = n u (n, A) + n d (fx, A) + n s (/i, A) , (67) 

while the pressure and energy density are evaluated through 



p(ji,A) = -b+ r dp 



\ djj, I d\i 



e(fi, A) = -P(/i, A) + /i [n u + n d + n s ) - /i e (/i) (n u - n e (//)) , (68) 

in analogy with the massless quark case, cf. Eqs. f)60p and (IBT]) . Note that the contribution 
coming from the electron number density in the above equations is rather small. 

In Fig. El we display our result for the total quark number density in cold deconfined 
quark matter as a function of the d quark chemical potential /i, and compare it to the cases 
of two and three massless quark flavors (for which beta equilibrium and charge neutrality 
require that \i e = for Nf = 3, but /i e ^ for Nf = 2). Upon comparison with the 
two-loop 0(a s ) result (not shown in the plot), we find that perturbation theory exhibits 
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FIG. 5: The total quark number density evaluated to 0(a 2 ) for locally charge neutral systems of 
2 and 3 massless quark flavors, as well as for the two light and one massive flavor case ('2+1'). 
All results are normalized to the density of three free massless flavors 3/i 3 /V 2 , and assume the 
values A^g- = 0.378 GeV, m(2GeV) = 0.1 GeV, while the renormalization scale takes the values 
3A/(/i u + fi d + Ms) = 1, 2, 4 (for N f = 2, 2A/( / u„ + fJ-d) = 1, 2, 4). As expected, the 2+1 flavor 
result matches the three flavor result at large \x and approaches the two flavor result at small \i. 



convergence for // > 1 GeV, in analogy with the massless case. Somewhat visible in Fig. [5] 
are kinks at the critical chemical potentials at which the strange quark density drops to 
zero, and below which the quark matter is net strange quark free. We suspect that this 
is simply a consequence of not having enough energy to produce strange quarks with a 
non-vanishing in-medium mass: The chemical potential is required to satisfy the condition 
ji > m mcdium (fi s ), where the parameter m medium (/i) can be evaluated by studying the effects 
of the finite chemical potential on the poles of a massive quark propagator. A study of 
the one-loop quark self energy at finite temperature was recently performed in Ref. jH3| (cf. 
[EH), and a simple generalization of these results to finite /i shows that for g^i -C M -C /i, 
the mass can be approximated by the formula 



m medium (A, fi) ~ ym 2 + -^M 2 + 0{a 2 s ) > m. (69) 

This leads us to argue that the chemical potential at which the strange quark density vanishes 
does not need to receive large non-perturbative corrections, as instead of confinement physics 
only energy conservation is involved in the mechanism. As a consequence, we can expect 
perturbative results to give quantitatively reasonable estimates for this critical chemical 
potential, at least if A > 1 GeV at this point (cf. the discussion in Section Hi Bl) . 

Moving up on the chemical potential axis, we note that for fi > 1 GeV, the strange quark 
mass becomes unimportant and one recovers the result of three massless flavors, discussed 
in Section HVBl Interestingly, studying the O(a s ) and 0(a s ) quark number densities, one 
observes the trend that the effects of the strange quark mass become less important with 
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increasing perturbative order. This indicates that accounting for the interactions accurately 
in fact makes the QCD EoS less sensitive to the quark masses, which we suspect may well 
generalize to further perturbative orders as well. 



V. STRANGE QUARK MATTER AND THE PHASE TRANSITION 

In this section, we will use the results derived above to investigate two mutually exclusive 
scenarios: The properties of absolutely stable strange quark matter and the onset of the 
confinement transition from deconfined quark matter to the hadronic phase. In both of 
these cases, we will evaluate the EoS for bulk matter (in contrast to systems of finite size), 
which is the relevant setting for the astrophysical applications considered in Section |Vl] In 
the strange quark matter calculation, our main goal will be to sweep the parameter space 
of the theory to find out if there is a region that allows for the existence of stable strange 
quark matter, while for the confinement phase transition case we will investigate if there are 
density windows that allow for a smooth matching of the quark matter and hadronic EoSs. 

For completeness - - and to better describe real world finite density quark matter in 
beta equilibrium — we will include in our EoS a contribution modeling the effects of color 
superconductivity (CSC). This is accomplished by adding to the pressure a term accounting 
for the condensation energy of Cooper pairs in the Color-Flavor-Locked (CFL) phase (see 



e.g. Refs. [22|, |59|, [60| ) , 

3vr 2 

where the baryon chemical potential is /i# = u n +Ud+Us and the gap parameter A approaches 
at asymptotically high densities the form [6lJ] 



Pcsc = , (70) 



A = - 6/ V e- 3 - 2 /^, (71) 

with b a constant. In this work, A itself will for simplicity be assumed to be a constant, 
which suffices to at least estimate the magnitude of the CSC effects. We will in each case 
study the values A = (corresponding to normal, unpaired quark matter) and A =100 
MeV, of which the latter can be viewed as an upper limit for the size of the gap. It should 
also be noted that while Eq. (ITUj) in principle obtains corrections due to the finiteness of 
m s , the condensation term has its largest value in the m s = limit. Thus, Eq. (ITUj) serves 
the purpose of indicating the maximal effect that quark pairing may have on the EoS of the 
system. 



A. Strange Quark Matter Hypothesis 

Stable strange quark matter configurations in vacuum can exist if there is a (strange) 
quark chemical potential \i for which the strange quark density in the system is non-zero, 
while the pressure is vanishing and the energy per baryon 

E/A= e(ji,K)/n B (ji,A), (72) 
is lower than for the most stable nucleus [l8[ ( 56 Fe and 62 Ni), 



E/A < 0.93 GeV. 



(73) 
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FIG. 6: Exclusion plots of stable strange quark matter, where the smallest possible value of E/A 
(required to be lower than 0.93 GeV for stability) is plotted as a function of the parameters of 
the theory. Left: Normal (unpaired) stable strange quark matter can only occur for comparatively 
large values of the renormalization scale. Right: The incorporation of the effects of the CSC gap 
significantly increases the allowed parameter space. Reducing the value of A^g has an effect similar 
to that of increasing A/ (fj, u + \x& + fj, s ). 



Here, the baryon density is related to the total quark number density n via ris = n/3. 
If stable strange quark matter were to exist, ordinary nuclear matter (made up of up and 
down quarks) would be only metastable, with a lifetime determined by the probability to 
generate strange quarks via several simultaneous weak interactions [62]. In particular, this 
would imply that nuclei with baryon number A > 6 would have a lifetime of ~ 10 60 years 



18[, consistent with current observations. In contrast, stable two flavor quark matter (zero 
strange quark density) is clearly ruled out by experiment, as it would imply an extremely 
small lifetime for ordinary nuclei because no strangeness needs to be created. 

In the following, we will consider the stability criteria for bulk strange quark matter 
(A ^> 1), which is both the relevant case for astrophysical applications as well as the most 
stable co nfig uration (finite size effects for small A only make strange quark matter more 
unstable [18|). We will also only consider quasi-static configurations, where beta equilibrium 
is maintained and the system is locally electrically neutral, of which the latter requirement 
follows from the fact that for strange quark matter no hadronic admixture can be present. 
Hence, we will require P( 2 '(fj,, A) + Peseta) = 0, where the three-loop pressure for 2+1 flavor 
quark matter, P@\ is taken from Section TlV CI and Pcsc from Eq. flTOj) . In addition, we will 
naturally enforce the condition that the density of strange quarks is non-zero. 

In practice, we sweep through the parameter space of the theory as follows. For both 
A = and A = 100 MeV, we vary the renormalization scale A, the MS scale A^g" as well as 
the strange quark mass m in the ranges 3A/ \fJ, u +^d+^s) = 1-.. 4, A^g- = 0.378+0.034 — 0.032 
GeV and m(2 GeV) ~ 0.100 ± 0.030 GeV. At each point of this parameter space, we first 
search for the smallest value of /i for which the non-zero strangeness condition n s (fi, A) > is 
fulfilled, which also corresponds to the minimal value of E/A. After determining the value of 
the latter quantity, we then choose B > so that the stability criterion P = is satisfied. It 
turns out that these constraints imply that our results do not probe physics below A ~ 0.95 
GeV, which gives us confidence in our method. We find that that for A = strange quark 
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matter is either unstable (E/A > 0.93 GeV) or violates experimental evidence (n s = 0) for 
most parts of the above uncertainty ranges, while a finite value for the gap parameter has 
the effect of relaxing the constraints significantly (cf. Refs. (63l . 64j for similar conclusions). 
These effects are summarized in the exclusion plots of Fig. 

The conclusions one can draw from the above results are clearly strongly dependent on 
whether one considers A = or A = 100 MeV. While the former case suggests a rather 
hostile parameter space for stable strange quark matter, perhaps offering an explanation for 
the absence of experimental evidence for it in direct searches [65N68|. in the CSC case our 
results are inconclusive. Thus, within the current uncertainties of the different parameter 
values, we can neither confirm nor rule out the existence of strange quark matter. Alterna- 
tively, if stable strange quark matter is indeed realized, then the restricted parameter space 
allows us to set limits on its properties, 4 as we will observe in the next section where we 
consider the masses and radii of dense stars made up of strange quark matter. 



B. Phase Transition from Quark Matter to Hadronic Matter 

If strange quark matter is not the true ground state of nuclear matter, then at some 
critical chemical potential (or range of chemical potentials) quarks must become confined into 
hadrons. In this section, we try to estimate the location of this confinement /deconfinement 
transition by matching our perturbative EoS for quark matter to existing results for the 
hadronic EoS. While it is clear that our framework can never capture the details of the 
confinement process itself, it is not excluded that our results could be used to get a reasonable 
estimate for the critical chemical potential. One should contrast the situation to that of high 
temperatures and small /i, where the matching between the hadronic and quark gluon plasma 
EoSs — the latter obtained through resummed perturbation theory — suggests a location 
for the deconfinement transition that agrees reasonably well with lattice data (see Fig. [7]). 
In the case of small T and high \i considered here, a quantitative test of our predictions will 
have to await further advances in non-perturbative solutions to QCD. 

To describe zero temperature hadronic matter in beta equilibrium above nuclear satura- 
tion densities, we consider three different EoSs, representing three classes of physical pictures. 
These are a nucleonic 'baseline' EoS by Akmal, Pandharipande and Ravenhall ((7l[, denoted 
'^4i8 + <5f + UIX*' in this reference), an EoS including kaon condensation by Glendenning and 
Schaffner-Bielich (j73|, denoted 'Uk = —140 MeV in this reference), as well as one including 
the effects of hyperons by Schulze et al. ([zl], denoted 'Vi8 + UIX + NSC89' in this reference). 
We have selected these three results because they constitute the most realistic and accurate 
calculations for the different physical scenarios, 5 in contrast to being 'maximally different' 



4 



It has been noted in the literature (see e.g. Ref. 28|) that for massless strange quark matter at high 



density, the EoS is extremely well (to the per cent level) approximated by the simple formula 

e = 4B cS + aP, (74) 

where B c ff is an effective bag constant. In our case, one sees that a similar ansatz works only in the limit 
of large energy densities, but breaks down where the effects of the strange quark mass (as well the running 
of a s ) become important. 

In the case of kaon condensation, one may argue that Uk — —100 MeV to —120 MeV would be more 
realistic values for the potential 0,0] > but we were unable to obtain a tabulated EoS for these cases. 
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FIG. 7: The matching of perturbative results for the quark matter EoS to hadronic EoSs at high 
temperature and density. Left: The pressure at \i = 0, obtained from resummed 0{o^J 2 ) pQCD [4| 
(using Eq. for a s with A^g = 0.378 GeV and Nf = 3) and compared with the result of summing 
up the effects of all hadron resonances with masses smaller than 2 GeV 32, 69]. For a cross-over 
transition, the pQCD results should match smoothly onto the hadronic EoS. The range in A for 
which this is possible corresponds quite well to the transition region determined by lattice QCD 
0, 70] (the grey band labeled 'LQCD' in the figure). Right: The T = quark matter pressure, 
taken from our O(aT) result with = 0.378 GeV, m(2GeV) = 0.1 GeV and B = 0, compared to 



three different hadronic EoSs 7l|, |72|] . The matching in general involves B ^ and does not have 
to be smooth, as the phase transition at zero temperature could be of first order. 



as in the selection criteria of Ref. [191 ] . 

The matching of hadronic and quark matter EoSs relies on imposing the following con- 
ditions: 

1. At the matching point, the pressure of the hadronic phase is equal to that of the quark 
matter phase. 

2. Both the hadronic and quark matter phases are locally charge neutral. 



3. The speed of sound c s = ydP/de has to be less than the speed of light in both phases. 

4. The energy density has to increase monotonically with baryon chemical potential 

The second of these criteria can _be easily relaxed by considering a two component system that 
is only globally charge neutral 
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This is, however, only a minor effect in comparison with 
the other uncertainties in the calculation, and we consider it misplaced accuracy to perform 
a detailed analysis of the two-component phase transition here. The last two criteria are 
on the other hand meant to impose naturalness on the resulting EoS: c s < 1 is required 
to maintain causality, and it would be quite bizarre if for any given /i matter composed 
of nuclei could have an energy density higher than that of quark matter (except if nuclei 
did not correspond to the true ground state of hadronic matter, cf. Section IV A)) . Together 
with the first condition, the monotonicity of the energy density also implies that above the 
critical chemical potential, the physical phase is always the one with larger pressure. This 
is an important consistency check of our procedure. 
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Case I ('low (J.') Matching Case II ('high p.') Matching 




M- B [GeV] H B [GeV] 

FIG. 8: The equation of state of hybrid quark matter, obtained by matching the quark matter 
and nucleonic EoSs at low (jib 0.39 fm -3 , Case I) and high (jib ^ 0.64 fm -3 , Case II) densities, 
respectively. The shaded regions are obtained by varying A, the integration constant B, and thus 
also the exact matching point, while keeping A-^ and m(2 GeV) at their central values 0.378 GeV 
and 0.1 GeV, respectively. 



In practice, we implement the matching as follows: For both A = and A = 100 MeV, 
and any of the three hadronic EoSs, we first fix A^g- and m(2GeV) within the ranges given in 
Section III and then vary A in its usual range to determine whether the criterion 1 above 
can be fulfilled for any (positive or negative) integration constant B and quark chemical 
potential. The successful cases are then subjected to the criteria 2-4, and the parameter 
values and hadronic EoSs that violate any of these are subsequently ruled out. As a first 
result, one observes that matching is only possible if 3A/(/i n + fid + Us) ~ 3.3. 

Summarizing the results of the matching calculation, we find that for both normal quark 
matter (A = 0) as well as for the CSC case (A = 100 MeV), matching to a hadronic EoS 
is only possible in two disjoint density windows: At low baryon densities, ub ^ 0.39 fm -3 
('Case F in the following), where all three hadronic EoSs are degenerate and above which 
the criterion 4 is violated, as well as at high baryon densities, % ^ 0.64 fm -3 ('Case IF), 
where matching is only possible to the purely nucleonic EoS and above which the criterion 
3 is violated in the hadronic sector. The matched EoSs of both cases are displayed in Figure 
[BJ from where we notice a significant decrease in the uncertainty of our results in comparison 
with the unmatched pure quark matter case of Figure [7J 

For quark chemical potentials above 0.4 GeV, the nucleonic EoS differs from the hyperon 
and kaon EoSs significantly. If one is to trust (one of) the latter two — and strange quark 
matter is not the true ground state of nuclear matter — our results suggest a confinement 
transition from quark to hadronic matter at or around the density of atomic nuclei ('case F, 
0.15 fm -3 < ub ^ 0.39 fm -3 ). 6 This is predominantly a consequence of our matching 
criterion 4 regarding the monotonic increase of the energy density, as well as the very high 



6 One should, however, note that the Case I matching can only be carried out in a rather small region of our 
parameter space, in particular requiring a large value for the renormalization scale A, and is thus in some 
sense less robust than Case II. We thank David Blaschke and Thomas Klahn for drawing our attention 
to this issue. 
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FIG. 9: The functional dependence between the pressure and energy density, given for a selection 
of different EoSs considered in this section. 



energies (corresponding to half of the Fermi pressure of a gas of free quarks) predicted 
by all of the hadronic EoSs at around % ~ 0.16 fm~ 3 . While keeping in mind that the 
matching process carries sizable quantitative uncertainties due to the perturbative nature 
of our calculation, we note that if the deconfinement transition at high temperature and 
low density — where e(T)/e^°\T) < 0.5 in the transition region |2j - - is any guide for the 
physics at zero temperature, our findings may after all not be so unreasonable. 



VI. ASTROPHYSICS APPLICATIONS 

Having obtained EoSs for both stable strange quark matter as well as hybrid 
quark/hadronic matter undergoing a phase transition, we are now finally able to compare 
our results to nature. At present, the only available observational window into the prop- 
erties of cold and dense nuclear matter are compact stars, which makes them the natural 
application of our results. 

Of special interest for us is the structure of a non-rotating compact star, in particular 
the relation between its mass and radius, because it is highly sensitive to the details of the 
underlying EoS of high density nuclear and/or quark matter. In addition, the mass-radius 
relation for a hydrostatic compact star is straightforward to determine in general relativity 
by solving the Tolman-Oppenheimer-Volkoff (TOV) equations (76| . 

dM(r) = 4nr 2 e(r)dr, 

= _G(P(r) + E (r))(Af(r)+4^P(r)) 
v ; r(r-2GM(r)) y ' 

Here, G = (1.22 x 10 19 ) 2 GeV~ 2 is Newton's constant in natural units and r the radial 
coordinate of the star, while the EoS enters through the function s(P) (cf. Fig. [9]). In 
practice, the TOV equations are solved by choosing a value for P at the center of the star 
(r = 0) and then integrating outwards to the surface, where the pressure vanishes. The 
resulting mass and radius can be calculated for any central pressure allowed by the EoS 
under consideration, and differ markedly between the different classes of EoSs. 

In Fig. [101 we show the mass-radius relations for several classes of compact stars, including 
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FIG. 10: The mass-radius relation for compact stars, obtained using A = (left) and A = 100 
MeV (right) in the quark matter EoS. We display the results for purely hadronic stars (containing 
only nucleons [7l| . nucleons with kaon condensation 73], or nucleons and hyperons 72]). pure 
quark matter stars ('strange stars', cf. Section [V Aj) and hybrid stars including both hadronic and 
quark matter (see text for details). Also shown in the plots are compact star mass observations 
from Refs. 
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sal- 



Purely nucleonic stars, obeying the EoS of Ref. [71 . 



Stars composed of hadronic matter including the effect from kaon condensation, 



cf. Ref. [73 



Stars composed of hadronic matter including nucleons and hyperons, cf. Ref. [72 . 



• Hybrid stars composed mostly of nucleonic matter and a small quark core, correspond- 
ing to Case II of the previous section. 

• Hybrid stars composed mostly of quark matter and a small nucleonic crust, corre- 
sponding to Case I of the previous section. 

• Strange stars composed of (stable) strange quark matter. 

These cases cover all the different realistic EoSs of hadronic and/or quark matter at small 
and moderate densities, and in particular include both the scenarios of stable strange quark 
matter and a confinement transition to the hadronic phase. Note that the branches for 
different EoSs in Fig. HO] could in principle be extended by allowing stars with higher central 
densities. These stars are, however, unstable with respect to radial oscillations [82|, |83| , so 
we do not display them in the figure 7 . 

For the Case I of hybrid stars, the curves shown in Fig. [10] corresponds to the maximal 
value of /i, for which matching is possible using the parameter values 3A/ + ^d + l^s) — 4, 



' The exception to this statement would be the so-called 'third family' of compact stars 8J, |85[ , which 
could occur if the matching between nuclear and quark matter takes place at densities above those reached 
in the center of the maximum mass neutron star. In our calculation, this possibility is not realized. 
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FIG. 11: Left: The effect of the variation of the baryon chemical potential /i pt , at which the 'case 
F matching between the hadronic and quark matter EoSs is performed (with 3A/(fi u + fj>d + Ms) 
I. = 0.346 GeV, m(2GeV) = 0.07 GeV). The value n pt = 0.99 GeV corresponds to the 
maximal value where the matching is possible. Lowering fi pt towards 0.93 GeV, the result smoothly 
approaches that of the strange stars, which are marginally stable with respect to neutron stars. 
Right: The effect of CSC obtained by varying the value of the gap parameter A. Relative to 
unpaired quark matter (A = 0), CSC generally decreases the radii of stable stars while the masses 
are hardly affected. 

Ams = 0.378 GeV and m(2GeV) = 0.10 GeV. Varying these numbers in the usual ranges has 
a large effect on the radius of the star, while the mass seems to be quite stable. The effect 
of varying the chemical potential at which the matching is performed is depicted in Fig. [TTJ 
For the Case II stars, the effects of the above parameter variations are considerably 
smaller, and our predictions thus more robust. The reason for this is that the structure 
of the star is most sensitive to the behavior of the EoS at relatively low densities, which 
in Case II is uniquely determined by the nucleonic EoS. One should, however, note that 
because we did not perform a detailed matching of the two-component mixture of quark 
matter and hadrons, including e.g. the effects of quark matter droplets in nuclear matter, 
our predictions for the radii of the resulting hybrid stars may contain sizable uncertainties. 
As discussed e.g. in Ref. [86[, the masses of hybrid stars we obtain should on the other hand 
be fairly independent on the details of the matching process, and can thus be considered 
accurate. 

Corresponding to stars containing only strange quark matter ('strange stars'), we show 
in Figure [10] the grey shaded regions obtained by varying A, A^g and m(2GeV) in the usual 
ranges. 8 The large size of these areas reflects the sizable uncertainty related to the behavior 
of the quark matter EoS at low densities. One should note that the maximal strange star 
masses are of the order of 2.75 M , with radii in excess of 17 km, in principle offering a 
clear signature for the detection of stable strange quark matter. The general effect of the 
color superconducting gap A is to somewhat decrease the largest possible radii of the stars, 
leaving the maximal mass basically unchanged. 

In addition to our predictions, we show in Fig. [10] observational data corresponding to 



Note that stars of this type have been considered in the literature before, using different lower order or 
resummed perturbative EoSs. See e.g. Refs. 33, 8^] and references therein. 
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compact star masses which allow one to exclude any EoS that does not pass through 

all the measurements. In principle, data on compact star radii could provide even more 
stringent constraints, but these are also much more difficult to obtain, and at present no 
radius measurement exists that is both precise and unambiguous. Within the next decade, 
there is nevertheless hope for the emergence of new data in particular for quiescent and 
isolated neutron stars 88] that could improve the situation significantly. 

From neutron star mass measurements alone, one can exclude two classes of hadronic 
EoSs, namely those of condensed kaons and hyperons without phase transitions to quark 
matter even at very high densities [72l.l73| (see Fig. fTOl) . Interestingly, our calculations for all 
types of quark matter EoSs (corresponding to both hybrid stars and strange stars) as well as 
for the purely nucleonic EoS lead to masses that are compatible with current observations. 
On one hand, this is unfortunate because it does not allow us to discriminate between the 
different physical scenarios. On the other hand, the fact that our treatment of the quark 
matter EoS naturally (and without any fine tuning) leads to compact star masses that are 
in agreement with observations is also rather pleasing. 

In light of the marked difference between our results (and hence conclusions) and those 



of Refs. [861 . l89j ]. some remarks are in order. The presence of a finite strange quark mass 
leads to a significant stiffening in the EoSs containing quark matter, which is not captured 
by either the usual MIT bag model or the ansatz of Ref . [86] , which can be seen to describe 
our EoS well only at \i > 1 GeV. For this reason, our calculation results in pure strange 
stars that can be more massive than pure hadronic stars, turning the argument in Ref. [89(] 
on its head (c/. also the discussion in Ref. 90]). Also, it seems to us that the largest neutron 



star masses displayed in Ref. [89| result from nuclear EoSs that can no longer be considered 
realistic. Therefore, we argue that a confirmed determination of a compact star with a mass 
in excess of 2 M & , as suggested in Ref. 89J, would seem to strongly favor — rather than 
rule out — the existence of quark matter in astrophysical objects. 

Clearly, our calculations can (and should) be improved in several ways, and also further 
astrophysical applications of our results are possible. Notably, a detailed calculation of the 
properties of the two-component system of quark and hadronic matter in hybrid stars can 
be performed with our quark matter EoS, which would allow us to quantitatively calculate 
also hybrid star radii. Furthermore, in this work the effects of quark pairing were modeled 
in a rather crude way, and one could certainly study them in a more detailed and realistic 
setting. Finally, our results can also be used as a basis for studying other astrophysical 
observables, such as a compact star's moment of inertia, as well as dynamical processes 



like compact star oscillations and cooling rates [9l|. We hope that our result for the quark 



matter EoS will lay the foundations for bringing the same rigor to compact star core physics 
than is today routinely applied to hadronic matter in their crust and surface. 



VII. CONCLUSIONS AND OUTLOOK 

Investigating the properties of cold and dense QCD matter with first principles calcula- 
tions is a notoriously difficult task. For densities relevant to real world physical systems, 
the strong coupling constant is not small, thus making perturbation theory converge at best 
slowly, while lattice QCD simulations are inapplicable due to the famous sign problem. Up 
to now, most studies of such systems have been performed using phenomenological models, 
most importantly the MIT bag model, where the effects of the QCD interactions are only 
visible in a single constant that is added to the grand potential of a system of free quarks 



33 



and gluons. If one aims at a quantitative description of the physics of cold quark matter, 
such a crude approach is clearly inadequate. 

In the paper at hand, we have tackled the challenge of performing first principles calcula- 
tions for zero temperature QCD matter using the machinery of perturbation theory. To gain 
better understanding of the thermodynamics of the system, and in particular of the effects 
of a finite strange quark mass, we have performed an 0(a 2 s ) evaluation of the equation of 
state of T = quark matter, keeping both the quark chemical potentials and the strange 
quark mass nonzero. We feel that such a calculation was long overdue for several reasons. 
One one hand, the existing 0(a s ) results available for the massive EoS exhibit a sizable 
dependence on the renormalization scale and offer very few possibilities for analyzing the 
convergence of the expansion — issues that can only be addressed by computing the next 
perturbative order for the quantity. In addition, in most of the practical applications of our 
results, the density of the matter is not large enough to warrant the approximation m s <C fi, 
thus leaving the door open for significant quark mass effects (33| . 

Our main result for the EoS of cold, deconfined quark matter with an arbitrary number 
of massless and one massive quark flavor can be found from Section IIIIEl and has been 
analyzed in Section IIVI As expected, our result for the 2+1 flavor case smoothly interpo- 
lates between the results of two and three massless flavors in the limits of low and high 
making the quark mass effects clearly visible. The ambiguity related to the choice of the 
renormalization scale is still sizable, but already at \i ~ 1 GeV clearly smaller than in the 
two-loop result. Investigating the effects of the strange quark mass, we find an interesting 
trend indicating that they become smaller with increasing perturbative order, suggesting 
that the EoS becomes less sensitive to m s when interactions are correctly accounted for. 

As physical applications of our result, we studied the (mutually exclusive) scenarios of 
stable strange quark matter as well as a confinement transition between quark matter and 
hadronic matter. Normal unpaired strange quark matter turned out to be stable only for 
large values of the renormalization scale A, in addition to which one needs to have relatively 
low values for both the strange quark mass and the MS scale parameter A^g. Incorporating 
the effects of quark pairing by adding to the pressure a term accounting for the condensation 
energy of the Cooper pairs, we on the other hand found that the constraints on stable strange 
quark matter were considerably relaxed. 

To study the hadronic-to-quark-matter deconfinement transition, we matched our quark 
matter EoS to various hadronic EoSs, including purely nucleonic matter as well as nu- 
clear matter in the presence of hyperons or kaon condensation. We found the physical 
requirements for a successful matching to restrict the behavior of the EoS considerably, 
discovering two distinct density windows where matching was possible: The 'low /i' region 
tib ^ 0.39 fm~ 3 and the 'high /i' region % > 0.64 fm~ 3 . Interestingly, we were unable to 
match the quark matter EoS to either the hyperon or kaon condensation EoS in the 'high 
/i' region. Thus, if we assume these hadronic EoSs to correctly describe nature, we must 
conclude that the deconfinement transition takes place at or around the density of nuclear 
matter n# ~ 0.16 fm -3 . 

Finally, we applied the various EoSs to describe non-rotating compact stars, investigating 
the mass-radius relations following from the different cases. Using reliable observations for 
the masses of the stars, we found that it is possible to rule out the cases of hyperon stars as 
well as stars with kaon condensation, in accordance with earlier findings Conversely, 



EoSs with a transition to deconfined quark matter can without any fine-tuning accommodate 
compact star masses close to 2M — close to the upper limit of confirmed observations - 
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and up to 2.75M for pure strange stars. Similar results for hybrid stars were also reported 
in Ref. [86 . 

While our evaluation of the 0{a 2 a ) EoS for zero temperature quark matter has evidently 
opened the door for a plethora of physics applications, it is clear in particular from the 
sizable renormalization scale dependence of our result — often overwhelming the uncertain- 
ties in the experimental values of m s and A^g — that an extension of the calculation to 
include four loop contributions would be desirable. Although the full al computation seems 
quite demanding, a meaningful and considerably more straightforward challenge would be 
to determine the coefficient of the 0(alln a s ) term in the expansion, which originates from 
the IR sensitive ring diagrams. This term would in fact only require a numerical evaluation 
of the two-loop gluon polarization tensor, a task of difficulty comparable to the calculation 
we have performed in this work. In addition to simply extending the region of validity of 
the current calculation, this result would allow one to apply new optimization schemes for 
the renormalization scale, thus hopefully significantly decreasing the uncertainty involved in 
its choice. 

Finally, the physics of cold quark matter can naturally also be tackled in various other 
ways apart from purely perturbative calculations. As we have tried to highlight in Section [TV] 
the calculation performed here can be used as an ingredient in several models, with an 
obvious example being the replacement of the simple MIT bag model by our three-loop EoS. 
To this end, we have attempted to make our results as accessible as possible, presenting them 
in terms of simple fitting functions above, and making all of them electronically available 
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Appendix A: The One-Loop Gluon Polarization Tensor 

Several parts of the derivation of the grand potential, presented in Section IHIt relied on 
various properties of the one-loop gluon polarization tensor at zero temperature and finite 
density. To this end, we will in this first Appendix review what is known about this quantity. 

The one-loop gluon polarization tensor is defined by the graphs of Fig. [3^. We divide the 
function into its vacuum (T — fi — 0) and matter (vacuum subtracted) parts through 

W»{K) = W v : c (K)+UZ t (K), (Al) 
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where we have suppressed the trivial color indices, noting that both parts of the tensor are 
proportional to 5 a b- A result that proves quite helpful in evaluating the VM graphs is the 
simple form, in which one can write the vacuum tensor, evaluated with iVj massless and one 
massive quark flavor. A straightforward computation leads to the result 
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Here, A and A refer to the constant and one- dimensional function appearing on the first 
and second lines of Eq. (IA3[) . respectively. 



For the matter part of the tensor, we on the other hand obtain e.g. from Ref. [49 



2tt 2 



,2 N f 

E 

i=i 



p 2 

dp — 



4e 2 - K 2 (K 2 + 2pk) 2 + Ae 2 p k 2 
1 H ^ — ; In 



arctan 



8p k 



(if 2 - 2pk) 2 + Ae 2 kl 



+ Aelk 2 - Ap 2 k 2 ^ 1 



(n mat )^(fc ,/c) 



9' 



2tt 2 ^ 
z/l i=i 



dp 



C P L 



8 e p p k ko 

2m 2 - K 2 , (K 2 + 2pk) 2 + Aelkl 1 
1 H — In p u 



8pA; 



(if 2 -2pA;) 2 + 4e 2 fc 2 j 



(A4) 



It. 



where the sum is over N[ massless flavors with ui = /Ji and one massive flavor with u n 
After a considerable amount of algebra, the remaining integrals in H^t may be evaluated 
analytically, giving (49 
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where $ = arctan^-, and in the sum it should be understood that m = for the N[ light 
flavors. It is also useful to consider the limits 
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which can be used show that the plasmon ring sum of Eq. fl53|) is both UV and IR finite. 

Finally, we note that upon separation of the contributions of the massless and massive 
quark flavors to II^t> one can decompose the combinations F and G, defined in Eq. fj48|) . 
in the forms 
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[cf. the Appendices of 



Appendix B: The Massive 2GI Graphs 

In this Appendix, we will go through the evaluation of the massive 2GI graphs of Sec- 
tion [TTIEnb] in detail, filling in the pieces missing from our earlier discussion. 



1. Scalarization 

The first step in evaluating the diagrams a-d of Fig. [2] is to perform all the Lorentz and 
color algebra associated with them and write them in terms of the scalar topologies T>i-T>^ 
of Fig. [12j The procedure can be automatized with the program FORM [93] , and leads to 
the results 
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FIG. 12: The topologies resulting from the scalarization procedure. Our notation for the functions 
T>{ from Eqs. (jBip onwards is such that a lower index at position k indicates the power of the 
propagator k in the figure, while the corresponding upper index gives the mass of the line in 
question, with c indicating that the line has been cut. For the uncut graphs, the chemical potential 
is present in each massive propagator, but not in the massless ones. 
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where the appearance of the #b m D a reminds us that it should eventually be replaced by 
the renormalized coupling constant. 

In the above expressions, we have used a compact notation, in which the functions 
Pi, P 2 , P3, P4 and P5 refer to the scalar topologies of Fig. [T2j In each case, a lower index 
at position k indicates the power of the propagator k in the figure, and the corresponding 
upper index gives the mass of the line in question, so that e.g. the two- loop graph of Fig. 
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with two massive and one massless lines, with the massless one raised to power two, would 
read 

4 I 1 1 O I = / d P0 / ,n W-T / d ?0 
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The chemical potential is naturally present in each massive propagator, but not in the 
massless ones, and therefore the massive scalar lines carrying \i will be referred to as 
'fermionic' below. Each scalar graph we encounter can be viewed as containing exactly 
one solid fermion loop, the orientation of which we are free to choose. We make the choice 
of following the numbering of the propagators in Fig. [121 so that A* i s always flowing in the 
direction of increasing propagator number. 



2. Cutting of the Graphs 

Next, we proceed to apply the cuts, introduced in Sec. IIII B b\ to the above 2GI dia- 
grams, and list the results obtained after performing one, two, and three cuts on the graphs. 
The cutting should be thought of just a bookkeeping tool, helping us to keep track of the 
various contributions to the grand potential. The correctness of the procedure can be verified 
by explicitly performing all the energy integrations in the graphs. 



a. Single- cut Diagrams 

Cutting always exactly one of the fermionic lines in Eqs. flBlj) - flB4j) . we obtain the fol- 
lowing results, where the symbol V indicates that the corresponding massive line in the 
diagram has been cut: 
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The terms with a mass derivative originate from integrals with squared lines, and in them 
it is assumed that m is set equal to m after performing the differentiation. In the last 
term of Eq. ( IB 8ft . the derivative is acting only on the integration measure of the external p 
integration. It should be recalled that in all the integrals appearing on the right hand side 
of these equations (including T> 5 ), the chemical potential has been set to zero. 

In deriving the above expressions, we have used the symmetry of the graphs in equating 
terms where a distinct, but topologically equivalent propagator has been cut. We have also 
used the FIRE algorithm [43] to simplify the final result by finding relations between the 
various two-point functions originating from the cuts. 



b. Double-cut Diagrams 

Continuing to use the same notation as above, we obtain for the double-cut diagrams 
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In terms where a squared propagator has been cut, it is understood that the momentum 
flowing along that line is P, so the external mass derivative on the fourth line of Eq. flB12j) 
only acts on the external p integration measure. 



c. Triple-cut Diagrams 

For three cuts, the two-loop graph in Fig. [2^ naturally vanishes, so we are only left with 
the diagrams of Fig. |2b-d. For these, we obtain 
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where it is again understood that whenever a squared propagator has been cut, the momen- 
tum flowing along that line is chosen as P. 



3. Collecting Results 



Simplifying the above formulae with the help of Eqs. ( 128]) and (I3T]) as well as Appendices 
[D] and [El we finally arrive at the following expressions for the one-, two-, and three-cut parts 
of the 2GI graphs: 
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In these expressions, the coupling g and the quark mass m now correspond to the physical, 
renormalized quantities, and the J n 's again denote integrals defined in Appendix iDl 



Appendix C: The VM Graph 



Here, we go through the evaluation of the vacuum-matter diagrams fl™ M and f2 VM in 
detail. Starting from the massive VM graph Qy M , we first perform the Lorenz algebra and 
obtain 
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where the tilde implies that the integrand in the scalar graph has been multiplied by the 
function fx prior to integration over P and Q. The four-momentum appearing in the argu- 
ment of f\ is naturally that of the bosonic propagator in the graph. 

The evaluation of the above integrals proceeds the same fashion as before, i.e. by applying 
zero to two cuts in the massive fermion lines. Dropping the \i independent zero-cut term, 
we obtain at the one-cut level 
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and for two cuts 

\m.2c 



n 



VM 



V 



-2d, 



(4tt) 



- (1-e) 



A 2e |2m 2 
fo 

2\e 







X(P-Q) 



((P _g)2 ) i +e 

m m 

c c e 



+ 



+ V. 



m m 

c c e 



(C3) 



To complete the evaluation of Eq. ( )32|) . we now finally plug in the values of the various 
integrals from Appendix [Ej arriving at 
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Next, we look at the cross term Oy M . Proceeding along the same lines as with 0™ M , we 
write 
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in which the first function is seen to vanish, while the second is given by 
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Here, the bar indicates multiplication of the integrand of T>^ by the function fi(m 2 /P 2 ) — 
/i(0), and just as before, the index c means that the line in question has been cut, the only 
difference being that now m = 0. 
In the end, one finds 
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An analytic evaluation of two of the remaining three integrals finally leads to the result 
displayed in Eq. ( )45l) . 



Appendix D: Integral Definitions and Results 

In this Appendix, we give the definitions of all of the integrals I n that have appeared 
in our discussion so far, as well as provide results for the ones needed in constructing the 
result in Section HTT| i.e. the ones not included in the functions Q n (m). When available, the 
results will be in the form of analytic expressions, while for the numerically evaluated ones 
we provide approximating pocket formulae. 
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Starting from the integrals needed in various parts of the 2GI computation, we have 
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where again a = (P — Q) 2 /m 2 and A(a) is as defined in Eq. (IC9I) . In these integrals, all 
four-momenta are massive and taken to be on shell (P 2 = —m 2 ), while /i only appears 
in the integration measures. When outside of the integral, the mass derivative again acts 
only on the mass appearing in the p integration measure (i.e. not on the ones in the q or r 
integration measures or on the masses inside the integrand), after which m is set equal to 
m. 

Finally, the integrals introduced in Section IIIIDI are defined by 
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Appendix E: Ad Integrals 



In this last Appendix, we complete our presentation by providing results for the various 
vacuum (T = u = 0) amplitudes that are needed as input in the single, double and triple 
cut integrals. They are given as power series expansions in e = (4 — d)/2 up to the order 
required, using the integration measure defined in the first section. 
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In simplifying the results, we have found the Mathematica package HypExp very useful |94 
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